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ABSTRACT 


An investigation is made of the drift motion of small particles 
under the influence of acoustic oscillations. The investigation is made 
to determine if the motion has a magnitude great enough to produce 
significant changes on the fuel distribution in the chamber of a 
liquid propellant rocket motor. The calculations are made for the 
motion in both a rectangular and a cylindrical chamber. 

In the rectangular chamber the gas oscillation is restricted 
to the fundamental transverse mode and motion in only one dimension 
is considered. The particle drift velocity, that is the non-oscillating, 
non-darnped term in the expression for the particle velocity is found 
in the solution of the second order equation. 

For particle motion in the cylindrical enclosure, only gas 
oscillations in the first transverse or sloshing mode is considered 
and motion ia restricted to a transverse plane of the cylinder. The 
particle drift velocity, again a second order term, is determined. 

The magnitude of the drift velocity is calculated using condi- 
tions found in a liquid propellant rocket combustion chamber. Dis- 
tances a typical fuel droplet would move during its average life time 
are calculated. The distances are small compared to the size of 


most rocket combustion chambers. 
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I. INTRODUCTION 

It is a well known fact that liquid propellant rocket combus- 
tion chambers often experience a phenomenon called unstable com- 
bustion. The gas properties, such as pressure, density and 
temperature, fluctuate during periods of unstable combustion. 

The fluctuations are not random and definite correlation exists 
between the fluctuations at different positions in space or different 
instants intime. Characteristically, the fluctuations have a fixed 
frequency and the amplitude of the variations can grow to values 
comparable to the steady state value of the property. 

The regularity and the large amplitude of the fluctuations 
distinguish the property variations observed during unstable 
combustion from the random variation in gas properties observed 
during "steady" operation of a combustion chamber. 

The observed frequencies of oscillation vary from 10 to 
as high as 10,000 cycles per second and depend on the driving 
mechanism or phenomenon and the chamber size. The frequen- 
cies can be divided into two general categeries: (1) a low fre- 
quency or chugging instability with frequencies in the neighborhood 
of 100 cycles per second and (2) a high frequency or screeching 
instability with frequencies in the range of 1000 to 10,000 cycles 
per second. | 

The low frequency phenomena are often caused by an 
interaction between the rates of propellant injection, combustion 


and rate of ejection of material from the chamber. 





ay 

In the high frequency case the variations of gas properties, 
during unstable combustion, are accompanied by motion of the gas 
in the chamber. In general, the motion of the gas and variations 
of properties are similar to the conditions which would be gener- 
ated by a strong acoustic oscillation. This observation suggests 
that the large amplitude variations observed during high frequency 
unstable Combustion are a result of a resonance which occurs 
when a driving mechanism, associated with the combustion phe- 
nomenon, is in phase with the natural acoustic modes of the com- 
bustion chamber. For example, in cylindrical chambers, the 
oscillations occur in the radial, tangential and longitudinal modes 
plus their various combinations. The first transverse or sloshing 
mode is the mode of oscillation most frequently found experimen- 
tally'*) (2) in combustion chambers with high frequency instability. 

The present investigation will be concerned with the high 
frequency instability. A complete and satisfactory answer to the 
question as to why this instability occurs has, as yet, not been 
found. In the attempt to answer the above question, many mech- 
anisms have been proposed that are based on the interaction of the 
pressure and temperature oscillations with the combustion process. 
The interaction depends on the existence of a time delay which is 
essentially the delay time in the response of the burning rate 
to a change in temperature or pressure. 

Since the delay times in combustion systems are very 


sensitive functions of the fuel-oxidizer mixture ratios, a systematic 
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variation of the mixture ratio in the chamber may set up conditions 
that are favorable for unstable burning. The suggestion has been 
made that variations in the fuel-oxidizer ratio can be produced by 
the relative mean motion of droplets of fuel and oxidizer under 
the influence of acoustic oscillations in the combustion chamber 
gas. 

This investigation is concerned with the initial require- 
ments for the production of the mean particle motion mentioned 
above. The requirements are first, that a particle drift motion 
does occur in a chamber with acoustic oscillations and, second, 
that the magnitude be large enough to result in the appreciable 
particle motion which would be necessary to produce significant 
changes in fuel composition in the chamber. 

The following two idealized problems will be solved. 
First, as a simple example, particle motion in a rectangular 
chamber with acoustic oscillations in the fundamental transverse 
mode is considered. In the second, more complex problem the 
calculation is carried out assuming a cylindrical enclosure with 
acoustic oscillations. The emall particle drift velocity induced 
by the first transverse mode of oscillation is determined. 

To calculate the particle drift velocity the equation for 
the velocity of the gas undergoing acoustic oscillations in the 
chamber is combined with Stokes drag law to calculate the forces 
acting on small particles. The resulting equation is linearized 
by introducing a perturbation expansion in powers of a term 


proportional to the amplitude of the drag force acting on the 





eas, 
particle as a result of the gas motion. It is necessary to go to 
second order terms in the solution for the particle position in 
order to obtain a mean drift velocity. 

Representative values of the drift velocities are calculated 
using conditions found in a liquid propellant rocket combustion 
chamber. Directions and distances the droplets would move, in 
their average life time, are determined. Finally, an upper limit 
is obtained for a maximum variation in particle density as a re- 


sult of the drift velocity. 
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Il. RECTANGULAR CHAMBER 

Development of the Equations. Before proceeding to the more 
complex and more realistic problem of the investigation of par- 
ticle motion in a cylindrical enclosure with acoustic oscillations 
it will be illuminating to look at particle motion in a rectangular 
enclosure containing acoustic waves. The rectangular case will 
illustrate the important features of the problem ina straight- 
forward manner since motion in just one dimension need be con- 
sidered. 

The method of solution for the particle drift velocity will 
be as outlined in the introduction. 

In solving the problem for the case of the rectangular 
chamber the following idealized assumptions will be made: 
1. The fluid medium is considered to be at rest and in a uni- 
form state except for weak isoenergetic fluctuations in velocity 
and state properties. 
2. The small particles are considered to be spherical and to 
experience Stokes law forces. 
3. The particles have negligible effect upon the gas motion in 
the chamber. 

Consider a rectangular chamber, as shown below, with 


acoustic oscillations in the fundamental transverse mode. 





ny = 2 
a. ee 
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Gas motion is restricted to oscillations in the mY direction. The 
linearized partial differential equation for the gas pressure dis- 
turbance in the chamber, under the conditions considered here, 


is a Simple wave equation of the form: 


3p ey 
yy? es ce Sh (2.4) 


The well-known solution to equation 2.1, restricting the solution 


to real values, is: 
ae 52) 
P = P sin 24 cos om (2. 2) 


Here, iS is the amplitude of pressure fluctuations of the gas in 
the chamber, C is the local velocity of sound, © is 27 times 
the frequency,” , and T isthetime. () must be Ve to sat- 
isfy the boundary condition that the gas velocity and, hence, the 
pressure gradient vanish at the chamber walls (Y suf ). For the 
first transverse mode considered here ™ is 1. The velocity and 
pressure fluctuations of the gas are related by the following equa- 
tion: 

e aU =— VAP - -9 (2. 3) 

ae yy 

The velocity of the gas in the chamber is given by integration 


of 2.3 as: 


UJ = fe Cos at Ge (aOE (2. 4) 





fe 
Equation 2.4 shows that the gas velocity is zero at the chamber 
walls and a maximum at the center of the chamber. The ampli- 
tude of the gas velocity oscillations is proportional to cos Te ; 
If the location of the particle is given by § , and if assump- 
tion Z holds, then the drag force on the particle is given by Stokes 


law and is: 


6m (U- Ft) | (2.5) 


The equation of motion for the particle is obtained by equating the 


inertial forces and drag force, and is: 


on SE, = Ciyhc (u- df ) | (2. 6) 


It is important to note that the gas velocity must be evaluated at 
the position of the particle, that is U = U(f,t) . If assumption 
3 holds, the gas speed,U , can be eliminated by substituting equa- 


tion 2.4 into 2.6; this gives: 


om JE. => ote Ie Ges re Sin wl ~ 6 Tu ob (2.7) 


This result can be simplified by the use of the following dimen- 


sionless variables: 
—— nN — ee << pt gta . 
Z wa a@ Z 3b | (2.8) 
Introducing the dimensionless variables in equation 2.7 gives: 


2 yo jae / : 
c= + k = = 8, C0S Hae sin 27 oO (2.9) 








Res 


Here, the constants are defined as: 


K= aamucd and = “Sémee tet (2. 40) 


The solution to equation 2.9 will give the position of the small 
particles, 2 , at any time, , while they are under the influ- 
ence of the oscillating gas in the chamber. 

Note that equation 2.9 is non-linear due to the appearance 
of Z inthe trigonometric function on the right hand side of the 
equation. Because of this nonlinearity, an approximate method 
of solution of the equation must be used. The approximation used 
depends on the fact that the term £& can be held to small values 
by limiting the magnitude of the viscosity, or other flow proper- 
ties appearing in Ss Note that this restriction on the calculation 
is independent of that required to obtain a linear acoustic equa- 
tion (2.4) in the first place. Thus, when the expansion in terms 
of € ‘ is continued to terms of higher order than the first, there 
is no obligation to extend the approximation in the acoustic equa- 
tion to higher powers of pressure or velocity amplitude. By re- 
stricting the analysis to arbitrarily small values of E , itis 
possible to linearize equation 2.9 by introducing a perturbation 


expansion of the form: 
ay) =a) che eee Oe (2.44) 


Substituting the expansion into equation 2.9 and separating the 
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equation into equations for the various orders of € gives: 








d“Z ‘dz — 2.42 

je + KG = 0 — 

fe P Kae = cos 42, SinatTd (2. 43) 
o / ; 

re + k d28 ae sin & Z, Sinawr Jd . (2. 414) 


In making the above substitution, cos 77 is expanded as 
e, / al e le 
cos ke = (cos¥e\I- 5 (2, 2+ EX xa] +---)- (Sina )EFAr€ Fz), 


For a particle at rest at J=0 and at the point Y , the 


appropriate initial conditions are: 


ne: dZ, a 
Z,(0o)= Zz , yy 69 = 0 





2 (0) « d2: (0) a (2. 45) 


\| 


Z,(0) = 2 (0) = O 


Solutions to equations 2.12, 2.13, and 2.14 contain terms 
which are constant, oscillatory, exponential and linear functions 
of the time. The purely oscillatory solutions to the above equa- 
tions are not of interest since no net transport of particles is 
involved and the solutions involving decaying exponentials have 
a negligible effect on the drift velocity after a few oscillations. 


Hence, the only parts of the solution that are of interest and that 
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result in a particle drift velocity are those solutions which are 
proportional to time. 

Looking at equations 2.12 and 2,13, it is seen that the so- 
lutions to these equations are of the form of constants, decaying 
exponentials, and sinusoidal terms. However, equation 2.14 
has Z,, multiplying the right hand side so that the product of the 
trigonometric solutions of Z, with [sinawo | will produce a 
constant term. The solution for 2. will then contain the first 
and presumably most important term proportional to time. 

In the following paragraphs the solutions for Z,, Z, 
and Zz, will be found. The chief aim of the analysis is to ob- 
tain the single term which is proportional to time so that the 
drift velocity, as defined here, can be evaluated. 

Solutions to the Equations. By inspection it is obvious that the 
solution to equation 2.12, with the given boundary conditions, 
is Z = Z = constant. 


Substituting Z,= Z in equation 2.43 leads to: 


+ kK 2 =E.C05 -— ae Sin Dipl (2. 16) 
; he 
The solution to the homogeneous equation is: 
ld Of 
Cae Orgs | (2.17) 


The trial particular solution is: 


e SinamJ +C, cos any (2.48) 








an 4s 
Substituting equation 2.18 into 2.16 gives: 


~ 47°C, sinarJ - 47°CcosarS + KarC,cosard Peo, 


—_k’ar, Sin 27 = jcos Ez sinanT |. 


Equating the coefficients of like terms gives the following values 


for the constants: 


cares C= Koos dz (2. 20) 
— 7 ) ~_ 7 ® 
Sie |S NERS ee 


The solution for = can now be written as: 


- Ror ihe. Z 
2,=[O.+q es SHE snare KSEE cosand], (aan 
— Tr = 


The initial conditions that Z, (0) = 21 (0) - 0 are used to eval- 
uate the constants C, and c . The following expression is 


obtained for the complete solution to equation 2.13: 


z= D7 Cos = Kk cos &? 2M COSze ee atk 
meat kK K2 7 +9 1° ~4m°k’— KF 
-}- cos ae Gin aes = bas Cosa Sh. 
-~40°-kK@ —k?2n-8 Wr (2.22) 


The differential equation for 2. can be found by sub- 


stituting 2.22 and Z, =Z into equation 2. 14: 


dt, Coes _ T cin TZ Dros Ez kes Zz . 
— + —_—— = ~~ am Sin z S45 a eee -f- > / = 
aa Kay ot 2” \4rkt k3 rs a sinamy 


Rly 


= 2160S ez po ; 
sin ¥2 (AT pOEE ) sinard 


m x sin EZ (55 cose | (sin ara) 


9. (2. 23) 


-F sin rr? Koos Fz Sumier dncos 27 | 


-K > 9r-3T, 








wo1Ze 
The above equation may be simplified slightly by using the follow- 


ing trigonometric identities: 
SIN oF cos$z7 = s+ Sin 72 
Sina2rJ cosarS = &sin4n7Tt (2. 24) 
(Sin ara) = SS - ~ COs ar J 


Equation 2.23 reduces to: 





/ G / : en og 
= z k 22s =\(F +G)sinawrd -Fe Gila yi 





+2 = ee CosfTrJ +£ S/n an | (2. 25) 
Here, 
(2 WSIn TE  ksinwe 
O * Sante? b= 16 (497+kK 9) ? 
po. Tisin 72 Gis kKsinre 
-2k (4074k ey” -& (494K @) 


The solution to the homogeneous equation is again: 
_kKY 
Ze = c. a Cs es 

The following trial solution will be used to find the particular solu- 


tion to equation 2. 25. 


es ={C,2 + C, Sin 7 Ob ae C, Cos ary 


2 Kae = asa (2. 26) 
+C,€°  sinawvT+ C,@ cosarT 


+ Cy Cos4ra + Cg sin 47 J | 





Loe 
Substituting 2.26 into 2.25, equating coefficients of like terms, 


and solving for the constants gives: 
7 ko WsSin wz CarpDekKeE a 
U+Ge Bk (41 renee) cd os. ey (ae 


B® | . 27° (Fre 
ote (ste) snare (22g) Sin 2rd 


2 5 -29rk 2 


KF 26) _aak =) | 
+ (Kee) COS ae Raa ee aye e aa 


a Z 
Sap erre ey “cos ar /. 


Here it is unnecessary to evaluate the constants C, and Gans 
since (C, is a constant and will disappear when the derivative 
with respect to time is taken to find the particle velocity and ‘oF 
is multiplying a decaying exponential which will not contribute 
to the steady state particle drift velocity. The non-oscillating, 
non-decaying term, a second order term, is Zz = Srv A: 
Using equation 2.11, the non-oscillating and non-decaying part 
of the expression for the particle position can be written as: 

ae pie OSAP ED ae (2.23) 
Or, using 2.8, equation 2.28 can be written in terms of the origi- 


nal variables as: 


es 2) (ome ies n & [224 eee 
Fo F + € (seas) (5 sE} E ) 7 (2.29) 
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Now taking the derivative of 2.29 with respect to time, the mean 


particle drift velocity is obtained: 


— me Pe Te | rf | 
ae oF mare "Oe aga ig ——_ 
Remember here that £ = Hemme fee and Ko Beet 


Equation 2.30 gives the steady state drift velocity of the 
particle which is at a point f in the chamber. The drift term 
is second order in og and depends sinusoidally on position. The 
greatest drift velocity occurs at an initial particle position of 
4 2tG. 

The fact that the drift velocity is second order can be jus- 
tified by the following argument. The drift velocity exists because 
of the gradient in the amplitude of the gas oscillations. Consider 
the forces acting on a particle due to the difference between par- 
ticle and gas velocity. The gas is oscillating at a fixed frequency 
but with an amplitude that varies spatially. Thus, during a 
single oscillation, the particle moves from a region in which the 
amplitude of gas motion is large to a region in which the amplitude 
of gas motion is 8™aller. During that part of an oscillation when 
the particle is near its maximum excursion in the direction of in- 
creasing amplitude of gas motion forces act to accelerate the 
particle toward the direction of decreasing arnplitude. Similarly, 
when the particle is on the opposite side of its mean position, 
forces act to accelerate the particle in the direction of increasing 


amplitude of gas motion. However, since the force acting on the 
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particle is proportional to the difference between the gas speed 
and the particle speed, the forces acting on the particle are higher 
in the regions of higher amplitude of gas motion. Hence, over a 
number of oscillations of the particle, a net force acts to accel- 
erate the particles toward the direction of decreasing amplitude 
of gas oscillations. The net driving force depends on the differ- 
ence between the forces acting on the particle on either side of 
the mean particle position and thus depends on the mean difference 
in the gas speed on either side of the particle mean position. 
Hence, the driving force and the resulting drift velocity is pro- 
portional to the product of the gradient in amplitude of gas motion 
and the amplitude of particle motion. Since both of these quan- 
tities are first order, the resulting particle drift speed is second 
order. 

The gradient in the amplitude of gas oscillation is pro- 
portional to sin? and the amplitude of the particle oscillation 
is proportional to cos 57 , 6g. see equation 2.22. Thus, the 
maximum drift speed would be expected at the value of Z which 
maximizes (cos Z7)(sinz 7) 7 = 1. 

Discussion of Solutions. It is of interest to determine the man- 
ner in which this drift depends upon the particle size. To do 

this, it is convenient to use the fact that e" can be written as a 

,_ -~4RK’ 


/ 
function of k , (¢ = Oe: 


nated from equation 2.30 to give the following expression for the 


Therefore, & “ can be elimi- 





drift velocity: 
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ax 
) 





CPS 
oe inthe , 
{V- OF 5 [gat S/ (2.357) 


To make the dependence upon particle radius explicit, write the 


eee eae 





particle mass as p7= CEN er” ; a ade Ce, re In these 
terms the drift velocity is given by: 
oe : Loe = 
4 (i 2 y 
RS ud Ses (2. 32) 


40°C? | (Beyer) + fs 


As the particle radius, o~, gets large it is evident from equation 
2.32 that the drift velocity approaches zero. If equation 2.32 is 


rewritten as: 


2 f Qeree 4 
‘ae fe ; ER = | ein Th (2. 33) 
a < C | (=f nes ft D 


it is evident that as the particle radius geta small, the drift ve- 
locity again approaches zero. Therefore, a maximum drift 
velocity exists for some finite value of particle radius. This op- 
timum radius, o. resulting in maximum drift velocity, can be 
found by examining the first and second derivatives with respect to 


co. The result of the calculation gives: 
& oO. rn. 
ae hic! ie (2. 34) 
Substituting the optimum particle radius, 2.34, into the 


drift velocity expression, 2.32, gives the following expression 


for the maximurn drift velocity: 
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> p* me 
(CO lef ) Sin 7 ‘ (2.35) 


In terms of the chamber pressure, Re , the maximum drift is: 


aes BY WE 
ap = sal sin | (2. 36) 


Note that while the optimum particle radius is a function of 





the chamber dimension, L , the maximum drift velocity, nr is in- 
dependent of the chamber size. Therefore, the drift motion is less 
important as the chamber size increases. 

Note also that the maximum drift velocity is independent of 
droplet properties and only depends on the gas properties through 


the speed of sound C , specific heat ratio h and pressure ratio 


BR 
re 
The ratio of drift velocity for arbitrary particle radius to 
the drift velocity for optimum particle radius is sirnply obtained 


from equations 2.32 and 2.34 and ia given by 


a alley | (2. 37) 
= [+ (Sen) | 


The drift speed drops rapidly for particle radius different from 
optimum radius. For example, for « = Zo or 26~* the ratio 


nw | 
Apx is about 5 
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Wy. CYLINDRICAL CHAMBER 

Development of the Equations. The techniques used in Section I 
are now applied to the solution of the problerm of particle motion 
in a cylindrical geometry. The problems are quite similar except 
that with the cylindrical geometry, (A) motion in two directions 
must be considered and (5B) centrifugal forces may be important. 

To investigate the motion of small particles in a cylindri- 
cal enclosure under the influence of acoustic oscillations, the 
same assumptions that were made for the rectangular enclosure 
will be used; namely: 
4. The fluid medium is considered to be at rest and in a uniform 
state except for weak isoenergetic fluctuations in velocity and 
State properties. 
2. The small particles are spherical and experience Stokes law 
forces. 
3. The particles have negligible effect upon the gas motion in 
the chamber. 

For the conditions given in assumption 41, it can be shown 
(3) that the gas motion is governed by the wave equation which is 


written below in cylindrical coordinates for the pressure distur- 


bance: 
/ be, = a Pa) 


Since we are interested in finding a particle drift velocity 
in a cross sectional plane of the cylinder, in analogy to the sirnple 


transverse waves considered in the rectangular chamber, only 
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oscillations independent of Z will be considered in this analysis. 
The solution to equation 3.1 for a circular cylinder with standing 


acoustic oscillations is given by: 


Pere 


P= PR Sin (ome) J. (4) | (3.2) 


ie 
Here, as in the previous problem, the frequency, )) , is ae , 


and, in this case the angular velocity, W. , is ee . The 
ae arise from the condition of zero radial velocity at the cyl- 
inder walls. The subscripts and ™ specify the modes of vibra- 
tion in the Y and “7 directions. At the walls “,,,, is a solution 


to 


a) ee a 


For each value of 77 there exists an infinite number of solutions 
to equation 3.3 since _J,, is an oscillatory function of its argument. 
Characteristic values of @%,,,, are listed in reference 3. 

For oscillations restricted to a transverse plane of the 
cylinder the simplest cases are the first radial mode (77=c ,/”=/ , 
m= 14.22), and the first transverse or sloshing mode (m=/,7=C , 

a 586). In investigating the possible particle drift velocities 
induced by the acoustic oscillations only the effect of the sloshing 
mode will be considered since it is this mode which is most often 
observed in liquid and solid propellant rocket systems experiencing 
high frequency instability. 

With the restrictions to oscillations independent of zZ and 


to real solutions, equation 3.2 can be written: 





Bris 


p cos¢? By (eae “| Gos Dimple (3.4) 


The coordinate systern used in this problem is shown in the follow- 


ing diagrar: 





In this case: 


VP = iy 8 (-siny) J, (73 “a se) cos arvt | 





+1 fr (cosp)(Cos 2rvt) (Te) a Gti) a) (nser)] (3.5) 


An expression for the gas velocity can be obtained from equation 
3.4 in a manner analogous to that used in the case of the rectangu- 


lar chamber. The result of the calculation is: 


os 


Y = ty aie Sigey oO) (ese Sin zrypt | 


Pe Re 
‘oh aearerrs cosp sinarvt |r (ig: - (2 aS )! (3.6) 


Fig. 2 shows the magnitude and direction of the gas velocity, 


as a function of Be , for various points in the cross section of 
Cc 


the cylinder. The calculations are made for a time such that 
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sin 27yt=1. Since the cross section is symmetrical, only the 
first two quadrants are shown in the plot. If can be seen in Fig. 
1 that there is a gas velocity node at ( equals 0° and 180° while 
the pressure node occurs at (P = 90° and 270°. 
Using assumption 2 the drag force on the particles is given 


by: 


6m (U-*) (3.7) 


—_ 


Here, as before, U is the gas speed and 7 is the particle speed. 
The particle position is given by is J. , the particle velocity by 


et = see g% ae lg hg a » and the particle acceleration by 
dt Gare Pe 


a= me f~ ot 
ga = 0, |e adesy] aif gdh. Be +AZe | 
jt oad 

The equation of motion for the particle can be obtained by equating 


the particle inertial forces to the drag force obtained from Stokes 


law. The result is: 


me St dey" i+t Ly om a2 Jt a Fh | 


Nel 
ery 





= 6 [z, 5 6 sin PT, (Tet *)sin arrl - p de | (oo) 


+b, a Fe ro Cos P Sin arvt/ Tr, (73 Teo?) - (ee) — 0 de) | 


In the above equation the same symbols are used for gas position 
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as for particle position since the gas velocity must be evaluated 
at the position of the particle (by assumption 3). Equation 3.8 
can be separated into two more useful equations corresponding to 
the a and Cg directions: 


Ln component 


dn 67uUe dre _ dy \* 
ee ee IE 2 (2) 


\] 


(3.9) 
“Sestak cos sin anyt [T,( Meat) (aia) 


= 
Lg component 


: ; TX, 
adh of +o Ly 4g =| BAER sine sin amyl T (mse | (3.40) 


Ags in the case of the rectangular chamber it will be conven- 
ient to nondimensionalize equations 3.9 and 3.10; the following 


dimensionless variables are used: 


ae oS Se D=. (3.44) 


Introducing the dimensionless variables into equation 3.9 and 3.10 


gives: 


ie component 
pe] 
z ae 
a, +k oS . z (42) 7 (3. 42) 


~ £7 %e cos sin amas’ | T, (rx, Z) - T, (17202) | 








— 


a Y component 


JZ Jo 
dT dT Aes 


E sing Sin 2TEA~ 2 5, (rot,2) | 


se 
(3.43) 


Here, €= sack and K = THe are parameters of the = 
lem and, except for constants, are the same as the > and K os 
the rectangular chamber problem. 

Equations 3.42 and 3.13 are non-linear due to the appear- 
ance of Z in the argument of the Bessel functions and are analo- 
gous to equation 2.9 in the case of the rectangular chamber. Since 

€ is proportional to BR » 4, and other flow properties, an argu- 
ment similar to that made about equation 2.9 is valid. The follow- 


ing perturbation expansion in powers of C is introduced to linearize 


equations 3.12 and 3.143. 
naa Gee eae & ae = (3. 44) 
oi) = (1) +E ¢,( + ORM +--- (3.15) 


In substituting the above expansions in equations 5.12 and 3.143, 

sin } is expanded as Sin ¢ = (sin $,)(I- £[9,+ EG] + ---) 

+ (cosh) (E ~, +€~$, - - --) and cos q is expanded as 

cos ¢ = (CoS fo)(I- E+E O,]°+---) - (singa)(Eg +64, - --) 

The Bessel functions are expanded in a Taylor series about the point 


Z, to get: 
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Bee a E ( WiyZo) - EZ, Woy T (TA Zo) -E°Z, rre,, J, (70, Z) 


10° 6 


252 772 
ZF (tro4,)" Te (rete%e) + SEER Ty, (1rete2) | 


T, (rev, 2) = IT, (tro6%) +ETx,2 T,(1x,Z) 7 a ea Gian) 


e532 
+E 1H, Z Te (14.2) — a T, (TTot,2,) + < (rota meee 7 ] J (7,2) 





E72" Trex 
— hoes ay (reiaes) | 





Ja (WeigZ) = 3 (Weg) FEZ Wor J, ( WeoZ)- a TL (we%%,) 


2 . o_o 
+&2, 1x, J, (1, %) - La T; (70%) + af (re,,) o(Wex,2) 





222 e,a 
om 3e Sree) al (1°, 2) 5 252 a fe (1, Z,) - 
re) o pe | 


The results of substituting 3.14 and 3.145 into equations 


3.42 and 3.13 and separating the results into orders of € are: 


_— 





ee component 
Pz. d& ~z/(db)* ~ 1 3.46 
yye +t KF z (5%) he 





JZ, Cian pane JG dh _ > (dbl = 
gy * Og rn a ay) ia 


Trate cig ATLAS cos Gy)=n), (np Z) + Ta (rai, ); | 
: 


a 
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QL 


J"% 


yaa t kgs -2nchae -2,(4f)- a7.38 db - 2, (ate) - 





eer 


Se SIN zrva J az TX%y_ COS P, J, (1Xyo ay 


— 


i == cos d, Ja (Win Z) + P Sin LT, (a5) - Bia al) (3. 18) 


l t component 


JZ. doe > J", : ar d& 
aoe I IF +2 oe PEO TS 





== 0) (3.49) 
IZ de 9h, dz d% | 2 7% Jb, 
27,483 tal ze * Mea + 22,558 Te ez ize raZz FR 


= Tr) SIND, Sin a | (3. 20) 


dm IPs Jd Pe dz, dt = dz J 
Ahk appli —— Ov 
ae, AZey Ar te Ty Me ee iy 1° 3 pe f 
-» fe a . a. _ wa 
“- ' i ge 7 J fe) / 2 - 
+a ek 1223, 78. +2e2, 08 1 acer i zs Mote 
Panto, obs x 2, te Se SIN P| Ltrx,5 y\ WZ » ene A (3,24) 
a ¢, COS Pf, 5 iva gk ) 510 Sra : 
If a particle is considered to be at rest at T=C andata 
point Z r p then the appropriate initial conditions are 
2.19) = Fy QO) = $, S8(0)= S% ia) = c | 
ZO) eS GY (0) = oP, (oO) = C | (3. 22) 
 —" ie -_ 
a) = ue = Gis = 2.3 ake ae : 
Se Ie 
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It is noted again that purely constant, damped or oscillatory 
solutions will not be of interest since no net transport of particles 
is involved. 

The form of equations 3.16 to 3.24 can be compared with 
equations 2.12 to 2.14. Since Q, , and Z, multiply the right side 
of the second order equations it is expected that any terms in the 
solution which are proportional to time will again appear in the sec- 
ond order equations 3.15 and 3. 214. 

The expectation that the drift velocity will be second order 
in € , specifically for the radial motion of the particles, also fol- 
lows from the physical arguments given on page 14. However, 
for the geometry considered here, centrifugal forces also act on 
the particle. Since the centrifugal forces are proportional to the 
square of the particle velocity and the particle velocity is first 
order in € , the drift which results from centrifugal forces is 
also a second order term, 

Solutions to the Equations. <A solution to the zero order equations 
3.46 and 3.49 is Z = constant = Z and d, = constant = o . The 
constant € is closely related to a drag force and can be called 

an interaction parameter so that to zero order in € there is no 
drag or interaction between the particles and the gas in the cham- 
ber with the result that the particle position is fixed. That is to 
zero order in € , the particles are not moved by the gas, but are 
at a fixed position. Using these simple zero order solutions, the 


solutions for Z,, Z, , @, , and @, can be found by the same 





e2le 
method used in solving equations 2.13 and 2.14 in the rectangular 
chamber problem. The complete solutions to equations 3.17, 3.18, 
3.20 and 3.24 are discussed in more detail in Appendix A. 
The first and presurnably most important non-oscillatory 
and non-damped terms appear in the solutions to the second order 


equations 3.48 and 3.24. These non-oscillatory and non-decaying 


solutions are; 





(Zn 
2 asc 4 (FE + DG) OF 
. K 


-amay 
(3.23) 
~-g(EH+DT) 


KE kt 2) J 


= a sin 


?, 


Here, 


oe ue LENS p | 5, (WXpZ) - Dy (1K,,; 2) | 
p= cos§ (nest, (ried) - Ete T(r9,2)] 
G = THesind | I, (TZ) - Ta (tre?) 

H = sin 6 | 1, BY (Tey) — er TT, (ro, 2) | 
a Cos ob PY, (Weouaas) 


Differentiation cf equation 3.23 with respect to time will 
give the particle drift velocity. The objection might be made that 
the decaying exponential terms in the solutions to the first order 


equations may contribute appreciably to the particle motion. A 
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calculation to determine the magnitude of the first order decaying 
exponential terms, Compared with the steady state drift velocity, 
is shown in Appendix B. In general, these terms appear to be neg- 
ligible for reasonable values of the parameters. 

Moreover, the terms which decay exponentially in time 
arise frorm certain assumed initial conditions of the particle posi- 
tion in a certain phase of the motion. Since there is no physicai 
way in which one initial condition or phase may be preferred over 
another, these can occur at random and hence the decaying terms 


can not contribute to the mean motion of large masses of particles. 
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IV. RESULTS AND DISCUSSION 
The Drift Velocity. Ignoring the oscillatory terms and decaying 
exponentiais the drift velocity of the particles can now be written 
down using the solutions listed in 3.23. Remember that the velocity 
of the particlesis A~ = 7 9% Ton DL. 


or in the dimensionless expanded form is: 


= T, c | die + ede + Ege | + Ty cl z, If 


Jv Jy 
(4. 4) 
aye 7 oY, doo do , pad , 
oy + €29h+ ER oe +6298 + oz, dh | 
The ca component of the drift velocity can be written as: 
2 pe 
= Ee | Bo ~ (FE +DG) or Be] BFE 8 ei 
& 2 =e : 
2 LK [kK + (2229*] K(k? + (202977 


in terms of the original variables. Note that (ap -- (FE Uap, 6) 
is a constant for any given 2 , and cylinder radius and does 
not depend on the particle dimensions. Let [& pr (FE+D 6) | =6 , 
then the “2 component of drift velocity is more conveniently 
written as 


E“c @ 
= 3 
“a = aki Ky (222a)] won 


a 


To see how the /2 component of particle velocity varies with par- 
ticle size the expressions for € and kk can be substituted into 


alg 


,, » this substitution gives: 
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Ci AR ee, a 
ee 
(ee) ee 


Vrm 


Or, since the mass of the particle equals SF pz aE o* , 4.3 can 


be written as 











4 Mh ee _ 
64D pariz 4g’ 
(Ue = oF ee 2 oe. (4.5) 


As o-->oco in equation 4.5,1%-—0 andas o->0O.U;>0O:; 
therefore, as in the case of the rectangular chamber, there is an 
optimum particle radius, oa”, which results ina maximum QV, . 
This optimum radius is again found by examining the first and 
second derivatives of the velocity with respect to o. The result 


of the calculation gives: 


Zo tet BOCs. 
on  «V4er ‘ane 


The same arguments can also be made forthe { component. 


Thus, the {Y component of particle velocity is 





2fis ; Se ee \ | 
ng = Gef-B leh +02) | oy e/—aeHe DT) | an 


8 [kfa?e (ARP) I} 2 fk lw PRY") | 


ee 
in terms of the original variables. Again, a EntDZ) is constant 
for agiven % , ¢ andcylinder radius. If A; (ei *DL) = es. 


then the { component of drift velocity can be written as 





a 
Since the expression for the /2 component of drift velocity (equa- 
tion 4. 3) differs from the (Y component of drift velocity (equation 
4.8) only by a multiplying constant, it is evident that the optimum 
radius for maximum Vo will again be ee Z6.9% ‘ 

Since the optimum particle radius is the same for the 2 
and ¢@ components of drift velocity, a particle of radius o”* will 
indeed experience a maximum resultant drift velocity. Note that, 
since = 2 in the case of the rectangular chamber, the optimum 
radius given in 4.6 agrees exactly with the optimum radius found 
in equation 2.34 for particle motion in a rectangular chamber. 

Again, in analogy with the rectangular chamber problem, 
it will be interesting to see whether or net the maximum particle 
drift velocity is independent of the cylinder radius, @ . The pare 


ticle drift velocity can be written as: 


~~, 


f pees . | 
"lnecna amar) “FS a arp] ? 


es a) 


Here @ and § contain the cylinder radius in the ratio ies and 
<7 

for a given ~< they are independent of the radius. It remains to 

be shown that © , kK and (2222) are independent of the cylinder 


radius. Since )) = ied , it is obvious that )) varies as the in- 


~~ 
me 


verse of the radius, and that en) does not vary with the 


fe 
radius. Since the optimum particle radius varies as a~ and the 
3 
mass of an optimum sized particle varies as q. , it is evident that 
e- S4eh ana p= 624¢% are both independent of the cylinder 
(Me YPCr (aga 


radius. Therefore, while the optimum particle radius is a function 


of the cylinder size, the magnitude of the maximum particle drift 


o3Z< 
velocity is not. This result is in agreement with the result ob- 
tained for the rectangular chamber problem. 
Another comparison with the rectangular chamber problem 
can now be made. If the optimum particle radius, equation 4.6, is 
substituted into the particle drift velocity expression, equation 4.9, 


the following is obtained: 


2 
es erp (EL Re Mes] te t0 

For (= On the constant § is © andthe particle motion is en- 

tirely in the radial direction. Fig. 1 shows that the gas velocity 

in the cylinder is also entirely in the radial direction when {% is om 
Since motion was restricted to one dimension in the rectangular geometry 
problem, a comparison of the maximum particle drift velocities is 


2) 
most meaningful if compared fora (@ of © in the cylindrical problem. 


When {= O is substituted into equation 4.10, the equation reduces to: 


es 6 fhe fk 
(W" = 5 ie ts seb eerie 


bo 


Se eR nS ay 
) "ax me (TX, aa) (iJ, (Toy 4 ~ 9 | M47) | AA. 41) 





For comparison, equation 2.36 is rewritten here: 


5 \a a 
ae x | | oan . (2. 36) 
The expressions are the same except for the terms on the right re- 
sulting from the geometry difference in the two problems. 

It is also noted that the ratio of particle drift velocity for 


arbitrary radius to the drift velocity for particles of optimum radius 


varies in exactly the same manner as given in equation 2. 37 for the 
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rectangular chamber problem. 
Magnitude of the Drift Velocity and Discussion. Thus far, with the 
assumptions listed in Section II] and the restrictions mentioned for 
€, the expressions for particle drift velocity are perfectly general. 

The expression for particle drift velocity, equation 4.9, will now 
be used to calculate the drift velocity of small liquid droplets ina 
liquid propellant rocket combustion chamber. Although the fluid 
medium in the combustion chamber is not at rest, the Mach number 
of the gas relative to the particle is usually low near the injector. 
Stokes law drag will hold if the droplet shape remains nearly spher- 
ical and if the Reynolds number of the flow is less than one which, 
for reasonable values of the parameters, means that the particle 
radius must be less than ten microns. The results of these drift 
velocity calculations will not be exact, but will give an indication 
of the mean motion the droplets experience in a combustion cham- 
ber. 

The following properties, for the droplets and gas in the 
chamber will be assumed: 
For the particle: 

@, = 2.4 slug/ft. > 

The particle radius will be picked to 

maximize the drift velocity and for 

the values listea here: 


Ms “5 
6 =2.75x1410~” ft. = 8.4 microns 


For the gas in the chamber: 
C = 2360 ft./sec. 


M2 4.54x 107° slug/sec. ft. 


& 3 


O 2A BAS AD * ato [te 
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y =4.3 
“T = 2500° K 
VY = 692 1/sec. 


For the chamber: 

QA=i it. 
In substituting the above values in the equations the form of © 
was changed to € = juss by substituting } a for oe 
Here, R and < are the mean chamber pressure and density. 

In calculating the magnitudes of the drift velocities, values of 
R of .4, .5, and .9 ft. were chosen for the initial particle posi- 
tion, and since the drift velocities are completely symmetrical, 
values of Y inthe first two quadrants of the cylinder cross sec- 
tion were considered. In particular, angles of ayer Zee 5° 45°, 
67. , 90°, and the corresponding angles in the second quadrant 
were used. Table I shows the results of the drift velocity calcu- 
lations. In both Table I and Fig. 2 the magnitude of the particle 
drift velocity is divided by ( Bey . Since the amplitude of 
pressure fluctuations in an unstable rocket chamber is a parameter 
which varies guite drestically, it is convenient to multiply the values 
listed in Table land Fig. 2 by the square of any desired ratio of 
pressure amplitude to mean charnber pressure to determine the 
actual drift velocity. 

Fig. 2 shows the directions and magnitudes of the drift 
velocities in the first two quadrants of the cylinder cross section. 
It is seen that the drift near the wall of the cylinder is greatest at 


(p = 90° which is to be expected since at this value of p there 
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is a pressure node and, as can be seen in Fig. 41, the gas velocity 
is entirely inthe ( direction. The drift velocities are generally 
in the radial direction, due to the centrifugal force exerted on the 
particles by the oscillating gas. The drift velocities are greatest 
when the initial particle position is at the mid-radius of the cylinder. 
Since the drift velocity is a second order phenomenon, it is difficult 
to find any correlation between the drift direction and magnitude, 
and the pressure and velocity distribution of the oscillating gas in 
the chamber except for that mentioned above. 

To obtain a range of values for the magnitude of the drift 
velocity, a value of .316 is chosen for 7a ; that is, a pressure 
amplitude amounting to 31.6 per cent of the mean chamber pres- 
sure. Then (Bey is 0.4 and with the values shown in Table 
I, the particle drift velocity varies from .39 to 4 feet per second. 

As mentioned earlier, the magnitude of the first order de- 
caying exponential terms in the solution of the cylindrical chamber 
problem are investigated in Appendix Bb. For this numerical 
example, the results show that after just two oscillations the mag- 
nitude of the velocity contribution of the deCaying exponential 
terms is leas than one-one hundredth of the magnitude of the second 
order drift velocity. 

It is also interesting to determine how far the particles 
move in the combustion chamber. In .01 seconds, which is close 
to an upper limit for liquid droplet life times in a rocket chamber 


(4), the droplets can move from .05 to .4% inches depending on their 
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initial position in the cylinder. This is a very small distance when 
compared with the size of most rocket chambers. 

In comparison, the drift velocity in the rectangular chamber 
will reach a maximum of i7.4 feet per second at an initial position 
in the center of the chamber. In .01 seconds the particles could 
move 2 inches which is considerably greater than in the cylindrical 
chamber, but still not significant when coniered with chamber size 
or wave length of the oscillation. 

Variation in Particle Density. Although it is not the aim of this 
present work to calculate in detail the variation in fuel-oxidizer 
ratio which results from the relative drifts of fuel and oxidizer par- 
ticles, it is of interest to calculate a few representative values for 
the change of particle density with time. The variation in particle 
density, 72 , can be computed directly from the continuity equation 
for nonsteady flows. The equation, in vector form, is 
In 4 7 (mH =O . (4.42) 
AG 
nr is the particle drift velocity andm=™7 aS h, (~) is the droplet density 
in number per unit volume. If the problem is restricted to the cal- 
culation of variations in density over short periods, then the varia- 


tion of 72 with time may be approximated by 


m : 
i UO), SIA Gea 19) oa 7 

Also, for short periods, the effects of evaporation of the particles 

may be neglected and hence the particle radius wili be constant. 


With the above restrictions equation 4.12 reduces to 
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re My (Gann) at : (4. 43) 


Values of wes are calculated for the chamber conditions given 

in the previous section. The calculation is carried out for the cyl- 
indrical chamber at As =Q.fand ¢ = 0 and for the rectangular 
Chamber at yp = 0.4. Atime of about 1/5 of the particle life tizne 
of 0.01 seconds is assumed and a radius equal to the optimum value 
for particle motion is used in the calculation. 

A 3 per cent reduction in particle concentration is obtained 
for the cylindrical chamber and a 10 per cent reduction for the 
rectangular chamber. These values of density variation are the 
maximum which can be obtained for the chamber conditions used. 
For the cylindrical chamber, the changes in particle density appear 
to be too small to cause a large change in oxidizer-fuel ratio and 
in fact the variations of the order of 5 per cent are much less than 


the variations which occur because of the injector design. 
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V. CONCLUSIONS 

It has been demonstrated that small particles in an oscil- 
lating gas field experience a steady state drift velocity in both 
rectangular and cylindrical enclosures. This result is based on 
the assumptions that (1) the fluid medium is at rest and in a uni- 
form state except for weak isoenergetic fluctuations in velocity and 
state properties, (2) the particles are spherical and experience 
Stokes law forces, and (3) the particles have negligible effect upon 
the fluid motion. 

The particie drift velocity is a maximum for certain opti- 
mum sized particles. The optimum particle radius is given by 
—— Van where ) is the frequency of acoustic oscillations 
in the combustion chamber. 

The maximum particle drift velocity is independent of the 
chamber size although the optimum particle radius is a function 
of the acoustic frequency and hence of the chamber dimension. 

The magnitude of the particle drift velocity in a typical 
liquid propellant rocket chamber varies from .4to 4 feet per sec- 
ond and in an average particle or droplet life time the actual 
movement will amount to just a few tenths of an inch. 

It is concluded that, for rocket chambers of reasonable 
size, the drift movement of particles has a small, if not negli- 
gible, effect on the distribution of fuel or oxidizer in rocket com- 


bustion chambers. 
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TABLE ] 
PARTICLE DRIFT VELOCITIES IN THE CROSS SECTION OF A 


CYLINDRICAL ENCLOSURE 


% 6 er ee GS 
(degrees) (ft/sec) (ft/sec) (ft/sec) 
ne 0 44.38 0 44.6 
ayy Fd ie: 10.6 -.3 10.6 
m| 45 i 5 -.4 7.6 
4 67.5 4.5 Sc 4.5 
4 90 - 3.9 0 3.9 
co 0 37.4 0 37.4 
a 259 34.4 -5.3 34. 8 
2 45 Stud -3.7 21.4 
a) 67.5 49.7 -2.6 19.9 
“fo 90 16.7 9 16.7 
9 0 10.9 0 10.9 
9 Le.) i2.2 -6.9 44.0 
29 45 1S. 4 “9.85 18.2 
39 67.5 13.25 26.9 19.8 
yy, 90 19-7 0 19.9 
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SOLUTIONS TO OUATIONS 33416 TO 3724 
In this appendix solutions are obtained for equations 3.16 


to 3.21 which are developed in Chapter III. 


= 
—, 


If a particle is considered at rest at J=o ata point Z, w 


then the initial conditions for the equations are: 


aS} er es Oe 
Z(o)= 2, $.(0) = , Ngee ee _ =, 
J 


48 ¢o O (4) 





QL 
VU 


2 (eye SE (0) = o, (Oy) = 





raz, UD) = 5250) = , (0) = a D> (p) = O 


Zero Order Equations. The zero order equations, 3.16 and 3.19 














arc. 

dZo JO, 2 4" ap IPo _ 
22.55 5G +R Ta + BOK TS : 7 
Z, 2% Bikar ces —- 4 > (dee) er.) (3) 


=, etconstant = Z and @, = constant = d satisfy equations 2 and 
3 and will be used as the solution. 


First Order Equations. The first order equations reduce to the fol- 


lowing upon substituting Z and b WitOos eu, amd 55.20: 


—_ 


‘y component: 


wet +7 Kae =|, (wa) Sin @ SiR arpar | (4) 
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(, COmponent: 





a , ES re = = _ 
PE eh gZs = (MY sin 2EPT Cos § [Tez] -T( maz) [| 


The solution for @, is obtained first and then the solution for Z, 


is found. 


Ee D 


a (i= eee 


WT 


(6) 


the Ly component equation may be rewritten as 


dd, dg, — , 2rva JS 
ya tT K ay Io SS seers 3 (7) 





The solution to 7 is determined by the usual technique of finding the 
complementary and a particular solution and adding the two. The 
constants appearing in these solutions are then eliminated by use 
of the initial conditions. 


The complementary solution to equation 7 is: 


Pt ao 
. = c + Cae (8) 
For the particular solution try: 
= 1 amypa Ss ampa J 
g = |Gz Sin 7 + C, Aer, (9) 


In this appendix, the subscript Cc will pertain to the comple- 
mentary solution of the equation and the subscript p will pertain to 
the particular solution. 


The first and second derivatives of equation 9 are: 
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Substituting 10 and 11 into equation 7 results in: 
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Equating coefficients of like terms gives: 























a 
—({aTVPa ~-k aT Vo = 
(ama) C, -K2t*= Cc, = D, (43) 
and 

27TVa _ (aTra \% = 14 

Katra c, -(2rPe\"C, = Oo. (14) 
The solution of equations 13 and 14 is: 
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The complete solution for q, is then: 
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The solution of 17 and 1418 is: 
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A similar process to that carried out above leads to the solution of 


Z, from equation 5. 
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then equation 5 reduces to: 
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Since this is exactly the form of equation 7 and since the initial condi- 


tions on p and Z are identical, the solution can be written down 


immediately by analogy and it is given by: 
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second Order Equations. The second order (, Component equation, 


3.18, reduces to: 
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when zero order solutions are substituted. 
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The values for ay , Z ,» and p, can be obtained from the solu- 


tions given in equations 20 and 23. Thus, from equation 20: 
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Substituting 29, 20 and 23 into equation 27 gives: 
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Equation 30 can be simplified by using the following trigonometric 

identities: 
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The results of substituting the above identities into equation 30 is: 
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Equation 32 is the linear second order differential equation which 
must be solved. The solution is obtained in the same manner as 
was used in the previous section to obtain values of the first order 
functions D and 7, . 


The solution of the complementary function is: 
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and the trial particular solution is: 
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The first and second derivatives of ee are: 
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Substituting 35 and 36 into equation 32, and equating coefficients of 


like term results in the following evaluation of constants: 
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Solving equations 37 to 44 for the constants gives: 
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The complete solution for rae can now be written as: 
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mo determine C and C, the initial conditions are substituted into 50: 
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Solving for C, andC, , relsults in: 
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The complete solution for Z, is obtained by substitution of equations 


53 and 54 into 50 and is: 
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The second order ly component equation can be treated 
in a Similar manner. Thus, substitution of zero order solutions in 


3.21 leads to: 
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Terms involving Z, and M, are evaluated by use of equations 20 and 


23 and are given by: 
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Substituting 62, 64 and 65 into equation 59 gives: 
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Using trigonometric identities, equation 66 is simplified to the 


following: 
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Substituting the above results into equation 67 and equating coeffi- 


cients of like terms we get: 
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Solving equations 72 to 79 gives: 
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The complete solution for ¢, is obtained by substituting the above 


constant into equation 69 and adding the complementary solution. 
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The constants C, and C, are found from 88 and 89 and are: 
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The solution for @, is then obtained from 87, 90 and 91 and is: 
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The terms proportional to time in equations 55 and 92 are 


circled and are the terms used to find the particle drift velocity. 
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APPENDIX B 
A DETERMINATION OF THE MAGNITUDE OF THE FIRST ORDER 
DECAYING EXPONENTIAL TERMS OF THE PARTICLE VELOCITY 
The first order decaying exponential terms in the solution 
for the particle position in the cylindrical enclosure are: 
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Here, 
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From equation 4.4 the expression for the particle velocity, contain- 
ing just the first order decaying terms, can be written as: 
a ce Sl + fg ce Ft | (2) 
dT] 
Substituting the expressions for Zz, and dq, into equation 2 the fol- 


lowing is obtained 
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For an initial position of the particle of 2 See eS, 4 = 45°, for 


oe = .316 and for gas and particle properties evaluated for the 
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conditions used in the calculations for Table I, the following ex- 
pression is obtained for the particle velocity due to the damped 


terins: 
— ; ‘z t 
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For T = .0023 seconds, which is time for approximately two os- 


; -~43500 _/0 ; 
cillations in the chamber, C =C ., and the velocity amounts 
to: 
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For these same conditions the steady state drift velocity as calcu- 


lated in Table I is: 
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Comparing equations 5 and 6, it is seen that after just two oscil- 
lations the velocity contribution of the decaying term is negligible 


compared to the drift velocity. 
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